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Abstract. 

In this paper the macroscopic Einstein and Maxwell equations for system, in which the 
electromagnetic interactions are dominating (for instance, the cosmological plasma before 
the moment of recombination), are derived. 

Ensemble averaging of the microscopic Einstein - Maxwell equations and the Liouville 
equations for the random functions leads to a closed system of macroscopic Einstein - Maxwell 
equations and kinetic equations for one-particle distribution functions. The macroscopic 
Einstein equations for a relativistic plasma differ from the classical Einstein equations in 
that their left-hand sides contain additional terms due to particle interaction. The terms 
are traceless tensors with zero divergence. An explicit covariant expression for these terms 
is given in the form of momentum-space integrals of expressions depending on one-particles 
distribution functions of the interacting particles of the medium. 

The additional terms in the left-hand side of the macroscopic Einstein equations for a 
relativistic plasma has murch in common with the additional terms in the left-hand side of 
the macroscopic Einstein equations for a system of self - gravitating particles (refer to. [1], 
[2]). 

The macroscopic Maxwell equations alsow differ from the classical macroscopic Maxwell 
equations in that their left-hand sides contain additional terms due to particle interaction 
as well the effects of general relativity. 
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1. Introduction. 

The idea of macroscopic gravity can be considered as an extension of Lorentz' idea (re- 
fer to. [3]), formulated first for electrodynamics, about the existence of two levels, micro- 
scopic and macroscopic, of understanding classical physical phenomena. Lorentz showed that 
Maxwell's electrodynamics is a macroscopic theory of electromagnetism, and the Maxwell 
equations could be derived from a system of microscopic field equations called now the 
Maxwell - Lorentz ones, by infinitisimal space - time regions averaging (refer to. [3], [4]). 

As it is known the macroscopic Maxwell equation for continuous media can be alsow 
obtained from the microscopic Maxwell equations by ensemble averaging the latter (refer 
to.[5]). 

The Einstein equations, whose right-hand side contain the energy-momentum tensor of 
matter, are phenomenological equations. It is natural to suppose that the Einstein equations 
(or their generalizations) for continuous media can also be obtained from the microscopic 
Einstein equations, i.e., Einstein equations whose right-hand sides contain the sum of the 
energy - momentum tensors of individual particles. However, due to the nonlinearity of the 
left-hand side of Einstein equations, the averaging of the microscopic Einstein equations is 
more complicated than one of the microscopic Maxwell equations (refer to. [6] - [10]). 

In (refer to.[l]) a method is developed for the ensemble averaging of the microscopic 
Einstein equations for interacting particles. ( We use the ensemble averaging procedure in- 
troduced by Klimontovich (refer to. [11, 12]) to derive the relativistic kinetic equation for 
a plasma. The same procedure was used by the present author in (refer to. [13, 14]) to 
derive a relativistic kinetic equation for a system of gravitationally and electromagnetical 
interacting particles in General Relativity accurate to within the second order for the inter- 
action constant.) This results to macroscopic Einstein equations for continuous media that 
are accurate to second-order terms in the interaction constant. The macroscopic Einstein 
equations for a system of interacting particles differ from the classical Einstein equations in 
that their left-hand sides contain additional components due to particle interaction. The 
components are expressed in terms of the two-particle correlation function of the particles. 

In (refer to. [2]) we got covariant expressions for additional components for the system 
of self-gravitating particles. The components are traceless tensors with zero divergence. The 
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expressions were obtained in the form of momentum-space integrals of expressions depend- 
ing on one-particle distribution function of the gravitationally interacting particles of the 
medium. The given expressions are proportional to the cube of the Einstein constant and 
the square of the particle number density. The latter relationship implies that interaction 
effects manifest themselves in systems of very high density (the Universe in the early stages 
of its evolution, dense objects close to gravitational collapse, etc.) 

The present paper is a direct continuation of earlier papers (refer to.[l], [2]), devoted to 
the derivation of the macroscopic Einstein equations for a system of self-gravitating particles 
to within terms of second order in the interaction constant. 

The objective of this paper is to obtain the macroscopic Einstein equations for system, 
in which electromagnetic interactions (for instance, cosmological plasma before a moment of 
recombination,) are dominating. 

The macroscopic Einstein equations for relativistic plasma differ from the classical Ein- 
stein equations in that their left-hand side contains additional terms due to particle interac- 
tion. The terms are traceless tensors with zero divergence. An explicit covariant expression 
for these terms is given in the form of momentum-space integrals of expressions depending 
on one-particles distribution functions of the interacting particles of the medium. 

The additional terms in the left-hand side of the macroscopic Einstein equations for a 
relativistic plasma has march in common with the additional terms in the left-hand side of 
the macroscopic Einstein equations for a system of self - gravitating particles (refer to. [1], 
[2])- 

The macroscopic Maxwell equations alsow differ from the classical macroscopic Maxwell 
equations for their left-hand sides contain additional terms due to particle interaction as well 
the effects of general relativity. 

2. Microscopic Einstein - McLxwell equations 

The method of deriving the macroscopic Einstein equations is discussed in (refer to.[l]). 
The notation we use here arc the same that in (refer to.[l]). 

Briefly, the method we used to obtain the macroscopic Einstein - Maxwell equations is 
the following. 

We start from the microscopic Einstein and Maxwell equations 



VfcF''^ = - — J\ (2) 
c 

Here G^^ is the Einstein tensor in a Riemannian space with metric Qij , x = Snk/c'^ 
is Einstein's constant (where k is the gravitational constant, c is the velocity of light). 
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T^^ is the microscopic energy-momentum tensor of particles, F^^ is the Maxwell's tensor, 
J* is the microscopic current vector of particles, T^^^^ is the energy - momentum tensor of 
electromagnetical field. Raising and lowering of indexes is accomplishment with the metric 
Qij , Vfc is a covariant derivative in a Riemannian space with metric cjij . 
The tensor T^^^^^ have the form 



Tii) = ^{ -F'lF'' + --g'^Fi^F^"^ ] ■ (3) 



The tensors T^^m) and J* has the form 



= E ^ac' I ^u\KNa{q\P.), (4) 



J' = T.^bC f^ulMq,P). (5) 
b V 9 

Here Ca is the charge of particles of species " a" , rria is their mass, g is the determinant 



of gij , p\ is the momentum of particles of spesies "a", = p\/ y gkjPaPa , 

d^p 

- is the invariant volume element in momentum spase [7]. 
Na{q^,PD - is the Klimontovich random function [11]: 

Na{q\P,) = E / dSSHq' - qli))S\P, - pf\§)). (6) 

Nere Ua is the number of particles belonging to species "a", s is the canonical parameter 
along the particle trajectories: ds = ^ gijdq'^dq^ ; g^^^ and p^''^ are the coordinates and 
momentum of the 1-th particle of spesies "a", which are found by solving the equations of 
motion: 

Ml) P\i) dp? 1 ~ ea~ 

-7^ = -p- = ^j,ikP(i)P{i) + —FikP^i)- 7 

ds rriaC ds rriaC y-' y c ^ ' 

Here Tj ik is the Christoffel symbol of the first kind given by the metric g^ . 
In view of Eqs. the random function (|^) obeys the equation 

P -K— + ^j,ikp'P -77Z- + —FikP -TTZ- = 0. 8 

ag* opi c opi 

Next we write the metric gij as 

9ij 9ij ~l~ hij^ (9) 
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and Fik as 

Fik = Fik + ^ik-, (10) 

Here gij = {gij) is the ensemble averege of the metric gij [1, 2], (Fik) is the ensemble 
average of Fik. Note that (hij) = and {uik) = 0. 

Parallel with the momenta p'^^i^ = macdq^^j^^/ ds we use the momenta measured in the 
metric Qij : 

P\i) = (^~\(1,p)pIi), a{q,P) = ds/ds = {gijpyY'^igikP^P^Y^'^ ■ (11) 

Here s is the canonical parametr introduced by g^j . 
The transformation from pi to Pi is given by 

Pj = QjkP^ = cygjkQ^'Pi- (12) 



THe Jacobian of transformation ([T^) , is (see [14]): 

«^^, (13) 



dpi I ^g 



dp 



where g is the determinant of gij . 

Now we introduce the function Na{q\pj) defined in the eight - dimensional phase space 
with coordinates {q,p) as 

Na{q,p) = E / dsS^q^ - q'^,){s))5\p, -pf{s)\ (14) 
1=1 •' 

where g^^^ p^j'^ are found by solving equations obtained from (|^) with the transformation 
(|T2|) taken into account (p* = g'^^Pj ). 

Note that the functions Na and Na are related in the folloing manner: 

Na{q,p) = ^,N,{q,p). (15) 
ga^ 



Equation for Na{q,p) can be obtained directly from (H) by replaicing the variables (^ 
and (|1|): 

pSt- + ^j,ikp'p^ + —Fikp" 



dq^ ' dpi c dp 
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d_ 

dpi 



fif.A^p'/ - -f^lj[,Aijp' ) Na . (16) 



Here 



k 

A,, = g,,-u,u,; = fi^- = f^- - T^- (17) 

/p^pi 



is the difference of the Christoffel symbols of second kind for the metrics gij and gij 

1 ^, „, 1 



a{q,p} ■ a{q,P) 



~lm rp (m rp 

y ^ mk y ^ mk' 



Im ] 



If in (0) and (H) we turn to the variables Pi and A^^ we get 



d^Pa 



a{q,p)J-uyiNa{q,Pa) 
g) ^9 



Ji f d^Pb [y i^r 



{q,Pb)- 



(19) 
(20) 



where d'^p/ y/—y is the invariant volume element in the unperturbed momentum spase. 
For subsequent calculation is it convinient to write the Einstein equations as 

p _i_ Y7 O™ X7 O"* -I- O"* O*^ O"* O" 



X^m^c^ / ^^a. /? (yikyjm - 7^9ij9km) u^^u'^Naiq^Pa) + xTi 



(21) 



Here Rij is the Richi tensor of the Riemannian spase with the metric yij , Vm is the 
covariant derivative in this spase. 

With the help (|^) and (|10D we can write the Maxwell equations (0) to within the first- 
order terms in hj 



+ V, {h!^UJ^^ - hluj^^) + \uJ^^Vkh = -AT^Y^^aj ^U\a^ (l - \h>i NM.Pa). (22) 



In (0), ^2|) and below when raising and lowering the indexes we use the averaged metric 
yij , h = h\. 

Let's expand the Einstein equations ( ^Tf ) up to the secand -order members in small hij 



and ujij : 



Rij + R^j^ + R[f H = Y^ X^aC^ J (yikyjm - ^yijykm ) ulu'^Naiq,Pa) + 



J \/ y 



(eZ)\(2) 



(23) 
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Here R^j^ is the sum of all terms of expansion Rij that are first-order in hij , Rlj' is the 
sum of all second-order terms in hij, ets.; L^j].^(/i) is the sum of all first - order terms in 
hij of expansion the expression 

[If- ~ _1~ ~ ^ 

T^j''^ is the energy-momentum tensor of averaged electromagnetic field Fij , (rj^j^^^ ^ is the 

sum of all terms of expansion Tj-f^ that are first-order in hij and cuij , {Tif''^^ ^ is the sum 
of all terms of expansion f^f that are second-order in hij and ujij . 
This expressions has the forms: 



(2) 



(24) 
(25) 



(l)m 



1 



g'"\-'^ihij + Wihij + Wjhu), 



(26) 
(27) 



-(1) 

-'ijkm 



1 /I 

"2 {hstu^u* + hst9^^) \ gik9jm - ^9ij9kr 
— (^-FiiFjmg''"^ + -gijFimF''^^ , 



rp(el) _ ^ 



PilPimh + ^hijFijjiF 29ij-^lm-Pkn9 ^ 



(e0\(2) 



(r.r*) 



47r 



u;iiFf-u;jiF,' + -F'^u;ir^ 



— FiiFj^nQst + -;j^9ijFikF^9st + -^9ijFimFst— 



' 2 -^.t Fkm9il9ji 



1 



FitgjWsm + FjiguQsm — 9ijFlm9st + -^gugjmFls iO ^h "*-|- 



(28) 
(29) 

(30) 



(31) 



2. Averaging of microscopic system of Einstein and Maxwell equations for 
the relativistic plasma 
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We average (^2]) and ( ]23|) over the paths (Ref. [11]- [14]) and introduce the one-particle 
distribution function 

Mq,p) = { [ ds5\q^ - <faii){s))6\p, - pf\s))) = ^{Na). (32) 
As a result we have the averaged Einstein equations in the form 

i?.,+A.,=x(T,f -(1/2)T(™)<?.,) + 

+ X (r^f - ^r(^o^.,) + +x (t!;^ - ^T«^.,) . (33) 



= Il^am^C^ / ^^^uluifaiq,Pa)- (34) 



Here 

is the macroscopic energy-momentum tensor of medium, T^'^j-) is the energy-momrntum ten- 
sor of macroscopic electromagnetic field (see. (^9]) , 

TJ?) = + ^^^^^'c^/^c.'™)- (35) 

is the macroscopic energy-momentum tensor of radiation in plasma, 

A., = - E X^ac' I ^{NaLi%^u-ut (36) 

« ^ a/ (7 



When obtained ( ^3] ) we taking into account that 

(4^) = 0, ((7t'y^^ = 0. 

Next we assum that 

Fikh'^j < Uij ets. (37) 

inside the correlation region. Consequantly ((^T/j'^y "*) is equal approximately to t/J-* . 
Taking into account the (^) - (28) we can write Ajj in the form 

Aij = Vk^'^.+fiij, (38) 



where 

2 



~{S'J■-S■K)P^^s, (39) 

= (/.rfi(^)'^s), (40) 

= {^'J', - s'K) Qhs + (41) 

QL = (4ffi.?'), (42) 
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+ -^g,,g^"' + u,uHf + u,u''5T - ^S^Sf} < Nahm > • (43) 

In (^) we reject the indices "a" on momentums and velosities of particles of species "a". 
The macroscopic Maxwell equations, obtained from after averaging, have the form 



VkF'' + Vkf"' + fi' = J\ (44) 

c 

where 

y.^^ = -(/^^O, (45) 

fi' = ^{u;'''Vkh) + X\ (46) 

X^ = -2nY.e, f -^ul,^{N,h), (47) 

'g 



74 

J' = Y.ebCn, [ -j^ulh- (48) 



b 

is the macroscopic current vector. 

To simplify still furthe, we only have to calculate hij and inside the region determined 
by the correlations radius and corresponding correlation time. Note, That distant collisions 
provide the main contribution to caculated macroscopic quantities. To consider this contri- 
bution it is enough to find hij and Uij from the Einstein-Maxwell equations linearised with 
respect to the metric Qij and the macroscopic electromagnetic field Fij . 

We assume the average gravitational field generated by the particles to be constant within 
the correlation region. In this case we can interpret gij within the correlation region as the 
Minkowski metric. We assume alsow, that the influence of macroscopic electromagnetical 
field on the microscopic field in correlation region is small (see (|37|)) . 

As a result we have linearized Einstein and Maxwell equations with respect the Minkowski 
metric . By employing the gauge VkY^ = 0, where jij = hij — {l/2)hgij, we get the 
foUoing form of the linearized Einstein and Maxwell equations 

□7^^' = - E ^Xm,c' f d^Mq, P'bWWl (49) 

b ■' 

= -ATI J2eJ d'p\Mq-P'b)< (50) 
b 

Here □ = g'^ViVj , ^b = Nb- n^h . 

Thus subsequent calculatios do not have a covariant form, but they are all done for the 
purpose of determining the components of the tensors ip^j , ^ij , (fij , /ij and t/J^ at some 
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(arbitrary) point (g) in the locally Lorentzian frame. In this reference frame the interval 
ds^ has the form 

ds" = dri" - {dq^f - {dq^f - {dq'^f. (51) 
The final result must be written in covariant form. 

The expressions for hij and we get from the Einstein equations (|i9|) linearized with 
respect to the Minkowski metric (which we still denote by Qij)) were found in ((refer to.[l]): 

x/ig)(r/,r/',p',k)$,(r7',q',p'), (52) 



^.^(r^,q) = E/ ^V/ d^q' j d^k J" di e-'^^"^-"^^ 



X 



xng)(r/,r/',p',k)$,(V,q',p')- (53) 

where (q = {q^,q'^,q^) is the three - dimensional radius vector in the given reference 
frame, and k = {ki, k2, k^) ), 

$6 = iVfe - Ubfb 

h^in, V',P', k) = --^{<u, - l9.,){e^'^^'--^ - e-'^^"'-)}, (54) 



-Ku''-lsi)kT]e-^'^('^'-^)}. (55) 



2 

In ( ^4|) and (^) the following vectors were introduced: 

kt = {k,k), k7 = {-k,k), 

where 



k = ^[{k,y + {k,)^ + {k,y] = \k\. 

Obviously, k^{k) = —k^i^—k). 

To obtain the additional terms Vk^^^ and /z* in macroscopic Maxwell equations we have 
to calculate the h = h\ and V^/i: 
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where 



where 



^ J J J J —OQ 

/iW(r7,V,p',k)<l>,(V,q',p'), (56) 
h^\r,,r,\p',\i) = !^f!^{e^Hv'-v) _e-^Hv'-v)\ (57) 
Vfc/i(r/,q) = E/ ^V/ d?(^ j d^k I" c/r/'e-^'^^'i-'i') 

b 

h.^ iv, V, P\ k)$6(r/', q, p) , (58) 
/^??(^, V,p',k) = f^,{kte^'^^'-^^ - k^e-^'^^'-^^}. (59) 



(27r)3fc 

Let's write the solution of (|50D in the form 



^^ik — diAk — dkAi. 

where [11]: 



Ai(r/,q) = E / A' / I d?k P dr]' e-'''^''-'^^ 

^ J J J J —00 

(//,//', p',k)$b(V,q',p')- (60) 

Here 

Af)(r/,V,p',k) = ^-!^«^{e^^(^'-'') -e-^^(^'-^)}. (61) 
The uOik have the form: 

^ifc(^,q) = E / A' / d^'q' I d'k r dr]'e-'''^^-^"^x 

^ J J J J —00 

Xa;2H^,^',P',k)$6(r7',q',p'), (62) 



where 



4^^, V,p',k) = ^^{(^.^4 - ktu[)e^'^^'-^^ - {k~ui - k,u[)e-^'^^'-^^}. (63) 
If substitute the (H), (|5|), (H), (HD, (H) to (||), (||), (H), (||) — we get the 

ij ' 



foUoing expressions for PJl , Ql^^ , Ajj , (/i^u;'"") , {uj'-'^Vkh) , Aj : 



= E / / ^^V / ^^'g' / t^'g" r dr]' r dr]" f d^k' f d^k" 

j^^J J J J J —00 J —00 J J 



X 
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Ql^^ = J2f d'p' f dY f d'q' I d'q" r dr]' f d^" f d'k' f d'k"x 

I J J J J J —OO J —OQ J J 



A.. = ^Y,^^^c^ j dy j dY J d^q' f dr]' J rf3fc'e-^'^'(^-^')x 



be 



+ u",u"'6r - \6^SJ'}hfl{v,v',p\^')nbnc9tcW;7^,q,p"), (66) 
(/i^,;'^™) = ^1 dV/ rfV'/ rf'q'/ rf'q"/' dv' f dr]" J d^k' J d'k" 



g-.k'(q-q')e-.k"(q-q");,^6)^^^ k')47(r7, k")n,n,^7,,(a;', x"), (67) 

(iv'^Vkh) = J2[ dY I dY' I d^ I d^q" r dr]' /' dr]" f d^k' [ d'k" 

, J J J J <J —OO J — OO J J 

DC 

g-,k'(q-q')g-.k"(q-q");,W ^, ^ J^,)^g k")n,n,^/,,(x', x") , (68) 

= _^27re, / / A" / ^^'q' dr]' f d'k'e-''''^''-'''^u"' 

, J J J J —OO J 

DC 

h^^^ iv, V', P', k')nbncgbcix'; r], q, p") . (69) 

In this expressions unprimed, primed and double-primed quantities refer to particles 
belonging to species "a", "b" and spesies "c" respectively. 

In (|4D - (^) we introduced the two -particle correlation function gab{x',x") (see Ref. 
[1], [2], [11], [12]): 

fabix, X') = fa{,x)fb{x') + gab{,X, x'). (70) 

Here fa{x), fab{x,x'), fabc{x,x',x") are the one-particle, two-particle and tree-particle dis- 
tribution functions respectively: 

( J ds5{x - Xai-S))) = fa{x), 

(/ ds6{x - Xa{s)) / ds'6{x' - Xb{s'))) = fab{x,x'), 



{J ds6{x - Xa{s)) J ds'6{x' - Xb{s')) J ds"6{x" - Xc{s"))) = fabc{x,x' ,x"), 
where 

S{x - Xa{s)) = 6\q' - q:{s))6\p, - p^{s)). 
11 



We denote the set of all variable (?7,q, Pi) by x, the set (V^q'jPD by x', while the 
momenta p'^^-^ are denoted by p' , and the p"^-, by p" . 

For the moments of random functions we have the formulas (Ref.[ll], [12]): 

{Na{x)) = njaix), (71) 

{Na{x)Nb{x')) = {naUb - naSab)fabix,x') + 

) 

+ riMaix) J ds'5{x' - Xa{s'/x)), (72) 

{Na{x)Nb{x')Nc{x")) = {nanbUc - UaUbSac - UaUb^bc- 
-naUcSab + 2naSabSbc)fabc{x, x' , x") + 
+ {nanc- na5ac)5abfac{x,x") J ds' S{x' - Xais' / x)) + 

+ {nanb - na5ab)5acfab{x,x') j ds" 5{x" - Xa{s" / x)) + 

+ {nanb - na6ab)5bcfabix,x') J ds"6{x" - Xbis"/x')) + 

+ nJab6bcfa{x) J ds'6{x' - Xa{s'/x)) J ds"6{x" - Xa{s"/x)). (73) 

Here Xa{s/x) stands for the particle path through point x of the phase space. Bearing in 
mind that = ^ i^afa and that fa is not a random function, we can easily obtain 
expressions for the avereges 

(iV,(x)<l>fc(x')), (iV,(x)$,(x')<fc(x")). 

In deriving (^) — (|69|) we assumed that ^ 1 and that x" = Xb{s"/x'), i.e. point 
x" is not on the path of particles of spesies "b" passing through the point x' of the phase 
spase. 

In work [2] two-particles correlation functions gab{x',x") are found for the system grav- 
itationally interacting particles. The two-particles correlation functions the for system of 
electromagnetically interacting particles were found by the author in Ref. [13], when getting 
the relativistic kinetic equation for the plasma (Eq. (18) from Ref. [13]). 

In our case we should find two-particles correlation function gab{x, x') caused by electro- 
magnetic and gravitationall interactions simultaneously. 

To obtain correlation function inside the correlation region we assume that /i^fc <^ 1 , 
cjife < 1 . Therefore fi^^- ~ fi^j^^ , ip^,, ~ u^^^ ■ 

After substituting (|53|), (|62D into (0), multiplying (pj]) by $b(x') and averaging we get 
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^{E/ dVlJ rf'k /_''^dr/"expHk(q-q')]: 



X 



(r^, V', k)pV"^A,, - ^...r (r^, r^', /, k)p'=A,,J (iV,(x)$,(x')$c(x")) 



(74) 



Next we assume that 



fabcix,x',x") ~ fa{x)fb{x')Jcix") 



In view of Eqs.(|7^), ( [73] ) the two partial correlation function gab{x, x') obeys the equation 

P'Tr-gab{x,x') + TijkP'p''T^gab{x,x') + —Fjkp''—gabix,x') = 



dpj 



dpj 



^{I "^^^'^I "^'"^"Z "^'^ f ^^"e^PHk(q - q")]x 

Min^ iv, V", P", k)p'p'" A,, - iv, v", P", k)/ 

c 

X fa{x)fb{x') j ds"5{x" - Xb{s" /x'))^. 



(75) 



In this equation for gab{x,x') we should put Tij^ = 0, Fj^ = 0, since we assume that 
within the correlation region the metric gij are constant and that the influence of macroscopic 
electromagnetic field is small. 

So we have the first order linear equation (|75|), whose right-hand side contains the sum 
of two terms. The first term caused by gravitationall interactions, the second one caused by 
electromagnetic interaction. 

Consequently, we can write the solution of Eqs. (|75|) in the form 



gab 



>, X') = 9ab \x, X') + 9ab \x, X'). 



(76) 



Here g'f,j^\x,x') caused by gravitational interaction, and (7^^''' (x, x') caused by electromag- 
netic interaction. 

The equation for g^^^\x,x') coinside with one in Ref. [1]. In [1] we got the gjf^^\x,x') 
in the form (see Eq. (46) in Ref. [1]): 



9ab {X,X) 



^ J-00 u 



xfi^i'V, r',p', k)ea;p[-^k(q - q') + ^(kv)(r7 - r) + ^(kv')(r' - r^')] + 



v' dr'r d 



00 p' 



dp'i 



%faix)x 



u 
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X nii:\T', T,p, k)exp[-zk(q - q') + -(kv')(V - r') + -(kv)(r - r/)], (77) 

c c 

where v = cUa/u^, v' = cu'f,/M'°, = {u^,u'^,u^), u'b = {u'^ , u''^ , u'^) . Here the 
subscript r indicates that after calculating the derivatives with respect p we must replace 
the arguments rj , and q by r and q+ -^(t — respectively. The subscript r' indicates 
that after calculating the derivatives with respect p' we must replace the arguments rj' , and 
q' by t' and q' + ^{t' — rj') , respectively. 

The equation for g^fl\x, x') coinside with one in Ref. [13]. In Ref. [13] we got the solution 
for c/lf (See (18) from Ref. [13]). 

With the preceding notation the result (18) from Ref. [13] takes the form: 

{el), n efeCe f (fk fV ,u'^ „^ dfb{x') r' 



oo 



x|44NP(-^k(q' - c(')sin{k{r' - r"))\ 



u'l d 

u'! d ~\ 
-^^NP(-^k(q' - c^')sin{k{T' - r"))] X 
u q^i J 

xexp[-{W){r]' - t') + -(kv")(r" - t]")]- 

x|44rb^P(-^k(q" - q').m(Mr" - r'))]- 
-44NP(-^k(q" - q').m(A;(r" - r'))] jx 



X exp\-{W){ri" - t") + -(kv')(r' - r]')]. (78) 
c c 

Here q^. = (r, g"") . 

After performing a differentiation with respect to ql in (|26|) we get a following expression 
for correlation function: 



-«>"=] (fc+e^'=(^"-^') - fc.-e-*'^^^"-^')) exp[-(kv')(77' - r') + -(kv")(r" - r^")]- 
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{k+e'''^^'-^'"^ - A;,-e-''=(^'-"")) exj9[-(kv")(V' - r") + -(kv')(r' - r]')]. (79) 



[-(kv")(V'-r") + -' 

c c 

It is evident, that the electromagnetic interactios in plasma are dominating. Consequently 



9ab {x,x) <^ g'^^'{x,x). 

That is why one can put the g^^^\x^x') instead of gab{x,x') in (|6^ - (|6^). 

If we now substitude (^) into ( p^ ) - ( |55D and integrate with respect q' , q" , k' , k" we 
get the following expressions: 



be 



X 



-A;, e*'^^^"-^') + kU-'^^^"-^'A exp[-(kv")(V' - r") + -(kv')(r' - r/')] + 

^ c c 



+ 



X 



-00 u"0 7-00 u'o ' ^ ' dp] 
[k+e^""^^"-^'^ - k^e'^^^"-^'^) exp[^(kV'){r]" - r") + ^(kv')(r' - V)]}x 

x/.r^^)(r^,r/',p',-k)(]^)(r^,V',p",k), (80) 



be 



"'"'ZC^^') I'^p' /^'^" /l'^' /l'^" I^^ 

m' dr' fr' dT''^„M{x'),,.„,^, ,,,,, 



X 



X (-A;, e*^(^"-^') + fc+e-*'^(^"-^')) exp[^(kv")(V' - r") + ^(kv')(r' - r/')] + 



X 



.-00 m"0 7-00 M'O'"' ' dp'] ' ^ ^ 

(fc+e-'^(""-"') - fcfce'^(^"-"')) exp[^(kv")(V' - r") + ^(kv')(T' - r/')]}x 

xfi^f\r/,r/',p',-k)n:^)(r^,r/",p",k), (81) 



6c 



■^UiUjg +-gijg +UiU dj + u^u d, - -d^dj \x 
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fd^k fv , , , J n' dr' r' dr" „.dfb{x') 

X {{u'u")6i - u'y) {k^e-''^^"-^'^ - kie'"^'"-^'^) x 
xexp[-{W){rj - t") + -(kv')(r' - ?]')]+ 

X (k+e-""^^"-^'^ - kre'"^^"-^'^) exp[^{kV'){r) - r") + '-{kv'){r' - 77')]}- 



X 



></:^''7:^''V^^{/:^/:^«^'')^i('''''"),-«Mx 

X (-A;,-e^'=(""-"') + A;+e-^'=(""-^')) exp[-{kv"){r]" - t") + -(kv')(T' - r)')]+ 

X (A;+e-^'=(^"-"') - kje""^^"-^'^) exp['-{kV'){ri" - r") + ^(kv')(T' - 7;')]}x 
— u"^ {{k'^e^'^^'^" — k"^e~^'^^'^"~^^^ ] 

^ L / x{ ^^^'(^- '^K" " *^' - """1 ^ 



X 



(-/t.-e''^^^"-^') + /t+e-^'^^^"-^')) ea;p[^(kv")(r/" - r") + ^(kv')(r' - 77')] + 
_ j^-^ikir"-r')^ exp['-{kv"){r]" - r") + '-{kV){r' - TyOjjx 



X 
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ixebelmbCnbUc 



be 



27r 



J d^p' J dY p dr]' J 



'^ f d^k 
1^ 



xu 



■ri) 



-ik(ri' —r])^ 



+ 



oo M'" J-oo W 

i 

— 

c 



u 



//O 



u 



/O 



^5) 



To simplify (^) — (^) still further, we proceed as follows. We assume that the distribu- 
tion function changes little inside the correlation region, so that in calculating the integrals 
in (pOD — ( p5D we can ignore, in first approximation, the temporal coordinate dependence on 
/ . We substitute the explicit expressions for hf^ and (Eqs. and into (|D|) - 
(|82D and evaluate the integrals with respect to r', r", 77', 77" and k. Then the expression for 
PJl becomes 



E 

be 



+ 



{u'u")5- - uy-]fb{x')^^Kfi{u',u") 



(86) 



Here we have introduced the notation K^j^^{u' ,u") and K^j^^{u',u") for tensors that in 
locally Lorentzian reference frame, in which gij = rjij is the Minkowski tensor, have the 
following form: 



k'S.(u',u") 



ffhf i,fl" ir'f dr"(e«l''-"- 

J rC J —CO J— 00 J —00 J —00 



xexp[^{kV'){ri" - t") + ^(kv')(r' - r]')], 



KflW. u") = ^ / ^ dr,' dr," dr" dr\e 
_^-ik{v' -v)^ gik{v" -v) _ e^^'^^^" ^^^) {k^ e~^^^'^" ^^''^ — k~e^^'''^"~'^'^)x 



'{Jk{v'-v). 
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xexp[-(kV'){r]" - t") + -(kv')(r - r]% 

c c 

After carrying out the integrals with respect r' , t" , rj' , and f]" these expressions take 
the form: 

Kf'{u,u ^— ' -^'''^ - 1-''^^ 



/d k ( 

^ {kc - kvyikc + kv")2 (/tc - kv")3(A:c + kv") j ~ ' ^' ^ ^ 

Kfliu', u") = -Kfliu', u") = -KfU^', u") (88) 

The above equahties hold only in a locally Lorentzian reference frame. To obtain covariant 
expressions for the tensors K'j^l^{u\u") and Kf^{u\u") , we take into account the following 
fact. The quantities K'j^{u\u") and Kpl{u',u") appeared in (|^) after the correlation 
function gab{x', x") was substituted to ( |5^ ) and result was integrated with respect to q', q", k' 
and k" . But the expression ( |79D for two - particle correlation function is a sum of two terms, 
which differ in that primed quantities referring to particles of species "a" are replaced by 
double - primed quantities referring to particles of species " 6" , and vice versa. It is after 
these terms were integrated with respect q',q",k' and k" that K^^l^{u' ,u") and K^^J^{u' ,u") 
appeared in (^61). Obviously, the both must be calculated in the same reference frame, for 
which it is convenient to take the center - of - mass reference frame, in which 

v' = V, v" = -V, = u'" = 1/^1 - vyc^ = 

. In this reference frame 



Here f = yff + fl + , where = t;" = u"" /u^ are spatial components of the vector v. 
A covariant generalization of (^) has the form 

The expressions for K^^^{u' ,u") and K^j^^{u',u") diverge as — > 0, i.e., for large impact 
parameters. The reason is that we integrate over an infinite region, while actually we should 
integrate only over the correlation region, where the metric is assumed to vary only weakly. 
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This difficulty is resolved, as well as in the case of kinetic equation deriving, by introducing 
a cutoff procedure in the divergent integral 

/■°° dk 
Jo F 

We set the lower integration limit to fcmm = ^/fmax , rather than zero, where r^ax is the 
size of the correlation region (the correlation radius). Then the above integral assumes the 
value l/2kli^ = (l/2)r^„^. 

As the experience of deriving the relativistic kinetic equation (refer to. [13], [14], [16], 
[20]) shows, more thorough investigations suggest that the integrals become convergent as 
r — » oo , with the contribution from the region where r > r^ax being infinitesimal. In Ref. 
[14], [16] there are estimates for Vmax in the case where the average metric Qij is the metric 
of isotropic cosmological model and in the case of gravitational interaction of particles. 

In the case of electromagnetical interaction of particles the parameter kmin is equal to 
^ , where r^i a radius of Debit, since electromagnetic interactions in the plasma are shielded 
under r > r^, . 

The tensor (0) possesses the following properties: 

K,, {u\ u") = Kij {u", u') ; Kiju" = Kiju'" = 0; K,, = Kj,. (91) 

Because of this the expression for P," simplifies considerably. The macroscopic Einstein 
equations incorporate not Pj", but the tensor ip'l- = —{l/2){6'^6j — SjSn)Pis- The expression 
for this tensor can be written as follows: 



^J,. = -(l/2)(5^5;-5j5:)i^: = E 



be 



-9) 



X^ebecmbmcnbricC^ r d^p' f d'^p" 



+ u'\u'u"){5lu'l + 5lu'') 



{u'u")Kjriu',u") U(x")^^ - fbix 



9) 



dp" 



Note that 



= 0, = 0, = 



(92) 



(93) 



Reasoning along similar lines, we can simphfy the expression for the tensor fiij , (pij , /i, 
which assumes the following form: 



be 



167r3 



x'^ebecmbmcnbUcC^ f d'^p' f d'^p" 



= [(^' + 7t)(«+ 
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X 



{z5f - Jr,m{n', «")/c(a:")^^, (94) 



be 



(95) 



Ilk n (I li\-C f ^ii\^ fb{-^ ) /nr\ 

u Jks{u,u )fc{x , 96 

opi 

Here z = {u'u") . 

In (p^, ( p6D we have introduced the notation Jrqm{u',u") for tensor that in locally 
Lorentzian reference frame have the form 

JimnW.u") = ^j^r di r d^' r' dr' r' dr^kte-^^^^'-^^- 

U'^U"^ J k-^ J -00 J~oo J-00 J-00 

-k+e-'''^^"~^">)exp\-(kV'){r]" - r") + -(kv')(r' - r]')]. 

c c 

After evaluating the integrals with respect 77', 77", r' and r" , we get 

T (1 ' 1 "\ — f ^ ^t^ti^n 

•Jlmn 5 ) 



u'^u"^ J A;3 (kv" - kv') I {kc + kv")^ 

h^h+h^ -I- h^h^h+ -I- h+ k+ -J- h+ -J- h+ 



{kc + kv")2(fcc - kv") {kc + kv")(A;c - kv")2 

^(fcc-kv")3/ ^""^ 

The symbol V.p. indicates that the integral is calculated as a principal value. 

Just as in the previous case, we specify (|97D in the center-of-mass reference frame, where 



v' = v, v" = -v, u"' = u"^ = l/^l-vyc^ = u'' 

. In this reference frame the components of Jimn{u' ,u") have the following form (the spatial 
indexes of three-dimensional velocity are lowered by using the tree-dimensional Kronecker 
symbol Sap) 

V f V V R \ 
JqOQ = -a{v) — , Jooa = -a{v)-^, Joa/3 = -a{v)SaP + I3{v) \ 5a(3 ^ j , (98) 
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J a 



Pi 



C 



Oaj r O/37 ^ :^ — 



+ 



07 



H ()/37 ^ 

C V J C 



(99) 



The function a and /? in (^) and ( P9| ) depend on the velocity v = ^vl + + v'^ only 
and have the explicit form 



TTC 



a 



2^ 1 + \ 



(1-$)' 



In 



1 + 7 



TTC 



2f (3 - 2g + 3$ 



+ 3 1 + 



In 



1 + 



(100) 



(101) 



Here we introduced the following notation for the integral 



1 



dk 



We set the lower integration limit to kmm = ^/fo ■ 

A covariant generalization of this results, which were obtained in the locally Lorentzian 
center-of-mass reference frame, to arbitrary reference frames has the form 

Jijk{u\ u") = A {giju'f, + giku'j + gjku'i) - z{gijul + giku" + gjku'-)- 



{u'iu''u'l + u'lu'm'l + u'lu''u'k) + ^zu'lu''ul 



, rf III I I I II , I II I , II I I \ , 



J 

II . I . I 



+ 



where z 



A 



kn 



[z-2) 



+ 



{2z - 1) 



C 



2nV2 

f^rmn 



(^- 1)2(2 + 1)1/2 (^ + 1)(2_ 1)5/2 

(z-6) (6z-l) 



In + 



+ 



(2-1)3(^+1)3/2 (2 + 1)2(2-1)7/ 



In (^z + \/ z"^ — 1 



The tensor Jijk{u',u") satisfies the identity 

Jijk{u',u")u"' -- 
Note that the tensor fiij is traceless: 

g'^Hij = 0. 



0. 



(102) 

(103) 
(104) 

(105) 

(106) 
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Let us now simplify the tensor t/J"* (see. (^51). Substitution of (|62D to (|35|) yields the 



expression for T^p . In view (^O]) - (|7|) we have (if ^ 1 ): 

= ^ (-^7.*^?,. + E / rfV/ ciV'/ rf^^g'/ rf=^/x 

X r dri' r dri" [ d^'k' [ rf3fc"e-'^'('i-^')e-^i^"(^-'i")x 

J —oo J —oo J J 

X uf \7^, r^',p', k')^(^)^'(77, V",P", k")nMc{x', x"), (107) 

Here we can not neglect by glfj\x',x") in expression (^) for ghc{x',x"). 
That is why 

t!;^ = rl^' + r^r^ . (108) 



One can get r^p if replace gbc{x',x") by {x',x") (see (|79D) in expression (|107|) and 
replace gbc{x',x") by g\^j\x' ,x") (see (|77D to obtain the expression for Tj^J^'' . 

Substitution (0) and (|7S|) to ( [TOTD yields the folloing expressions for t^^^ ^ and p^ ^ ' 



^" 4- 2f2vr)% J .fT-r^J .f7Z7P ' ' ^ ^''^ + 



-(0 

6c 



2(2vr)%7 yp^y yp^' 
+ ((5/ + 6ju"i)u'^ — g^'^ {u'iu" j + M"j'u'j)]x 
X « - <..-) J(J(«', «")/c(x")^^- (109) 



[gr) V- Xf^becfribmcnbricC^ f dY f dY' ro ^Prg 



pq_ 



-(5>", + Sy,)u'' + g^\u\u", + n'>',)]4^;V, u")f,{x") x 

^^-.{/,(x')[ (z^ - ^) + + ^) - 2zu\u'f]}. (110) 



x 



Note that the tensors r/f^-* and r/j-* are traceless: 



^..^'; = 0, ^7^^7^^ = (111) 
In ( [T09[) and (^) the tensors 4ppu',u") and ^^(^^^(m', m") have the form (^), where 



A and i? have the forms (|103|) and ( p.04| ) respectively. But in the expression for J^pqiu'^u'' 
we must put kmin = ^/rn, where td is the radius of Debit, since the electromagnetic 
interaction in plasma are shielded under r > td . In the expression for Jl,l'-^{u\u") we must 
put kmin = ^I'^g^ where is the radius of correlation for gravitational interaction. As the 
experience of deriving the relativistic kinetic equation (refer to. [13], [14], [16], [20]) shows, 
more thorough investigations suggest than the integrals become convergent as r — oo , with 
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the contribution from the region where r > Vg being infinitesimaL In Ref. [14], [16] there 
are estimates for in the case where the average metric Qij is the metric of isotropic 
cosmological model and in the case of gravitational interaction of particles. 

4. Macroscopic system of Einstein and Maxwell equations for relativistic 
plasma 

As a result were obtained the macroscopic Einstein and Maxwell equations in relativistic 
plasma. They have the forms: 

G., + + /X., - xr\f = xT^j, (112) 

471 

VfcF'^ + Vfcvp'^ + /i' = (113) 

c 

Here Gij is the Einstein's tensor of the Riamannian space with macroscopic metric Qij , 
F^^ is the macroscopic tensor of electromagnetic field (Maxwell's tensor), J* is macroscopic 
current vector, T^- is the macroscopic energy-momentum tensor. The last is the sum of 

(m) 



macroscopic energy- momentum tensor of medium T^j iuS)^ energy momentum tensor of 



macroscopic electromagnetical field T^f^ (B^) and macroscopic energy-momentum tensor 



T^j^ (109) of electromagnetical radiation in plasma. (In cosmological plasma in the last case 



one should say about the energy - momentum thensor of relict radiation.) 

The Einstein equations of the gravitational field for continium media, obtained here, 
differ from the classical Einstein equations by the presence of additional terms VfcV^fj , fJ^ij 
and —X'^ij'^'^ in the left-hand side. It caused by particle interaction. The forms of this 



tensors are (p2D, (Q) and ( |110| ). The third term, {—X'^ij'^^) is the addition to macroscopic 



energy - momentum tensor of electromagnetic radiation, caused by gravitational interaction 
which multiplying on x ^ind moving from the right-hand side of macroscopic equations to 
the left-hand side. 

The macroscopic Maxwell equations differ from the classical Maxwell equations by the 
presence of additional terms Vk<^^^ + /i* • The additional terms in Maxwell equations caused 
by particle interaction and by effects of general relativity. 

The tensors V^v^f^ , fiij , t^j"^^ , V^v^'^* and /i* are expressed in (|92D — ( pBf ) and ( |110|) 
in terms of one-particle distribution function specified in the eight-dimensional phase 
spase in which all four components of momentum are independent. The transition to the 
seven- dimensional distribution function Fa{q\pa) is made according to the formula 



nafa{q\Pj) = Fa{q\Pa)S{\/g^'^PiPm " m^c) . (114) 

Here the function Fa depends on the spatial components of momentum only. Greek indexes 
are used to denote spartial components. 
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By integrating (|92D — (|96D and (|110|) with respect to p'g and p"q we can whrite down 
the tensors Vk^'^,, Hij, t^-''\ Vk^^^ and /i* as 



9 9 9 ^ 



be 



+ u"'{u'u"){6ju': + 6{u';)\ {u'u")Kf^{u',u") (f,{x")^^^ - F,{x"^^''^''" 



dp'a 



be 



9 9 9 ^ 



d^p' f d^p" 



dp'o 



(115) 



ik V- Xebelrribm^c^ 



be 



271 



d^p' f d^p 



p"J{-9) ' P"y{-9) 



dp'a 



{u'u")Kfa{u\u")x 



(116) 



dp'a 



dp'a 



(117) 



be 



XelecTnbmlc^ 



4tt 



d^p' f d^p" 



(gr) ^ Xebeemlmlc^ 

'''' "t- levr^ 



u"'jUu',u")F,ix"f^'^'''^ 
d^p' 



[{{u'u")6i-u'y^)]x 



dp'a 
d^p" 



(118) 



p"^i-g) ' P"yi 



-{5tu", + 5jn",)n'^ + g^\u',u", + «",«',)] J^^7(m', n")F,(x")x 
d 



dp'n 



y.{Fb{x')[ [z' - ^) 5i + (;.^ + ^) - 2^..\«'^]}. (119) 



Here 



and 



are the invariant volume elements in tree - dimensional momentum spase of particles spesies 
"b" and "c" respectively. 

The greek index a in ( p.l5| ) - ( |118| ) takes the values 1,2 and 3 only (the spartial index). 
The derivative with respect to p'„ in ( |119| ) should by calculated as all four components 
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of momentum are independent. The dependence of p'^ on p'^ is taken into account after 
differentiation with respect is completed only. 



In (|TT8D and (|TT|) the tensors and 4f}{u' 



u") have the form ( p.02| ) , where 
A and B have the forms (|103| ) and ( |104j) respectively. But in expression for J^p^(u',u") 
we must put k^in = l/rn, where vd is the radius of Debit, since the electromagnetic 
interaction in plasma are shielded under r > rn . In the expression for 4l'^{u\u") we must 
put krain = ^/f'g, whcrc Tg IS the radius of correlation for gravitational interaction. As the 
experience of deriving the relativistic kinetic equation (refer to. [13], [14], [16], [20]) shows, 
more thorough investigations suggest than the integrals become convergent as r — » oo , with 
the contribution from the region where r > Vg being infinitesimal. In Ref. [14], [16] there 
are estimates for Vg in the case where the average metric Qij is the metric of isotropic 
cosmological model and in the case of gravitational interaction of particles. 
The tensors (p^j , /i^ , t^j"^^ /i* must obey the additional conditions 



0, 



Vifi' = 0, 



(120) 



:i2ii 



since the divergenses of G 
Equations (|12C 



VkF , VfcV? , J' vanish. 



121) impose some restrictions on the parametres r^) and Vg dependence 



on the coordinates and the relative velosity of particles. The latter can be expressed via of 

z 



iu'u") 



The macroscopic energy-momentum tensor t/^™^ of medium and the current vector J* 
can by alsow written in terms of one-particle distribution function as follows: 



S'p 



pi. 



-9) 



PiPjFaip), 



a p^ 



d?p 



P'FM- 



122) 



:i23) 



The system of equations ( |112| ). ( [L13| ) must by augmented by the kinetic equation for 
in relativistic plasma. In the case when the electromagnetical interaction og particles are 
dominating, the equation for Fa was derived in Refs. [11], [13]. 

The covariant form of this kinetic equations is 

i^/a , -p dfa e^j fc^/a 

+ ^j,ik^ + — <Fik>U — = 

ag* opi c opi 



E- 



(124) 
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where 



X { - gij[{u, U'f - 1] - UiUj - u[Uj + {u, u'){UiU'j + u'iUj)} 



(125) 



with {u, u') = u'-u^ . Primed and nonprimcd quantities refer to particles belonging to 
species a, and b respectively and L is the Coulomb logarithm [15] 



5. Conclusion 

The macroscopic equations of the gravitational field in relativistic plasma differ from the 
classical Einstein equations by the presence of additional terms 



on the left-hand side due to partial interaction. 

These terms are proportional to the square of Einstein constant and to the square of 
particle namber density. 

The macroscopic equations of the electromagnetical field in relativistic plasma differ from 
classical Maxwell equations by the presence of additional terms 



on the left-hand side due to particle interaction and due to effects of general relativity. 

This terms are proportional to the first power of Einstein constant and to the square of 
the particle number density. 

Hence these terms can play an important role in continuous media of very high density 
only. Such density are possible in the early stages of the evolution of the Universe and inside 
objects hear gravitational collapse. Therefore, it is natural to look for applications of the 
derived equations primarily in the theory of early stages of the Universe evolution and in 
gravitational collapse theory. 




(126) 



Zij = Vk^-j + liij - xr-j 
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